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Abstract 

We compute non-perturbatively the renormalization coefficients of scalar and 
pseudoscalar operators, local vector and axial currents, conserved vector and 
axial currents, and O^^^^ over a wide range of energy scales using a scal- 
ing technique that connects the results of simulations at different values of 
coupling p. We use the domain wall fermion formulation in the quenched ap- 
proximation at a series of three values of /?, 6.0, 6.45, and 7.05, corresponding 
to lattice spacing scaling by factors of two. 

I. INTRODUCTION 

Lattice QCD has proven to be a powerful approach to calculating from the first principles 
the mass spectrum, weak decay constants, weak matrix elements, quark masses, and hadronic 
structure functions. Such studies require calculating matrix elements of composite operators. 
In order to extract physical predictions in the continuum limit from Monte-Carlo simulations 
in lattice QCD, one in general needs to determine the normalization of these operators in 
terms of physical requirements. 

In calculations involving the operator product expansion one needs to know the Wilson 
coefficients. The corresponding calculations are done perturbatively, and therefore must be 
carried out at a high enough energy scale. To evolve those coefficients down to the scale at 
which lattice calculations are performed, one usually also uses perturbative methods. 

A possible approach to renormalization of operators defined on the lattice is lattice per- 
turbation theory. This approach, though, suffers from many difficulties, the main being poor 
convergence of the perturbation series. One source of this difficulty was identified by Parisi 
and Toe [||], and later by Lepage and Mackenzie 0, and is due to the presence of tadpole 
diagrams. The authors of [Q] proposed a method to improve the convergence properties 
of these series which is known as tadpole-improved perturbation theory. Nevertheless, the 
calculations using this method are complex, and therefore are rarely carried out beyond the 
one-loop level. Since decreases only logarithmically with decreasing a, the systematic im- 
provement of perturbation theory is very difficult, leaving the possibility for large systematic 
uncertainty in extraction of physical results. 
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In some limited cases renormalization of lattice operators can be obtained from Monte- 
Carlo simulations using the cliiral Ward Identities. The examples include vector and axial 
vector currents, and the ratio of pseudoscalar and scalar densities. Unfortunately, this 
method cannot be applied to general composite operators for which there is no corresponding 
Ward Identity. 

A method for full non-perturbative renormalization of lattice operators by means of 
Monte-Carlo simulations was proposed in The procedure mimics the approach taken 
in continuum perturbation theory. The renormalization conditions are imposed directly on 
quark and gluon Green's functions, in a fixed gauge, with specified, external lines carrying far 
off-shell momenta. This defines the so-called Regularization Independent scheme (RI). At 
high momenta the results can be directly related to the perturbative calculations, provided 
that corresponding conversions are made in order to match different renormalization schemes 
used in continuum calculations (usually the MS scheme). 

A necessary condition for the RI method to work is to have the off-shell momentum scale 
of the external lines, /i, much smaller than the inverse lattice spacing and much larger than 
Aqcd- The former condition restricts the value of reliable momenta from above for a given 
lattice spacing a. The lower bound on momentum comes from the need to compare with 
perturbation theory which is only applicable at an energy scale above the non-perturbative 
regime, 3> ^qcd- These two constraints leave a pretty narrow range of momenta that can 
be used to calculate renormalization coefficients. 

To extend the range of momenta over which an operator is renormalized, one can perform 
simulations at different values of coupling constant j3 = with overlapping regions of 
reliable physical momenta. By combining the results of these calculations, it is in principle 
possible to relate renormalization coefficients of a given operator between quite different 
momentum scales fi. By starting with sufficiently high momenta where the perturbation 
theory can be trusted, operators can be renormalized at a series of scales down to a few GeV 
where most lattice calculations are done. 

To make the matching between results obtained at different /5's meaningful, it is impor- 
tant to insure that it is done under identical physical conditions. In particular, mass and 
volume dependence have to be taken into account. 

The difficulty that arises from the finite volume effects can be approached in several 
ways. A direct way of taking the infinite volume limit turns out to be inefficient due to 
the strong volume dependence at low momenta. In this work we instead keep the volume 
explicitly finite. The matching between calculations at different /3's is done at the same 
physical momentum and physical volume. This procedure is described in detail in Section 

In our simulations we use the domain wall fermion formulation 0. This approach 
is computationally more demanding than the usual Wilson fermions, but offers significant 
advantages that are more than enough to compensate for the additional computer time 
required for the calculations. The DWF formulation provides an easy procedure of taking 
the chiral limit, independently of the continuum limit. No additional fine-tuning is required. 
By working at sufficiently large lattice extent in the fifth direction, the residual quark mass 
arising from the chiral symmetry breaking coupling between the walls can be made very 
small. Then the limit of zero bare quark mass corresponds to the chiral limit with high 
accuracy. We use this property to eliminate the mass dependence of our results, by taking 
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the chiral limit. 

The excellent chiral properties of the DWF formulation suppresses the mixing of such 
operators as 0^£^^ with operators of opposite chiralities, Ref. [Q. This is a serious compli- 
cation for studying the renormalization of such operators using fermion actions with large 
chiral symmetry violations. 

It is important that in the DWF formulation the chiral properties are improved not only 
for the on-shell calculations, but also for the off-shell ones. This is one of the reasons why the 
DWF formulation is ideal for using the non-perturbative renormalization method described 
in the next section that imposes a renormalization condition on the matrix elements of 
operators between external off-shell quark states. Another important consequence of these 
good chiral properties of the lattice theory with the DWF is the absence of order 0(a) errors 
in the results for the matrix elements of operators. 

We use the quenched version of the DWF method in our paper. For an extensive analysis 
of the properties of this formulation the reader is referred to the early works on this subject, 
Refs. |] and 0. 



The paper is organized as follows: In Section ^ we explain the method of non- 
perturbative renormalization. The scaling technique is described in Section |T|. In Section 
|V] we present the results of numerical calculations. 



II. THE NON PERTURBATIVE RENORMALIZATION METHOD 

In our calculation of non-perturbative renormalization coefficients we follow closely the 
procedure first proposed by the Rome-Southampton group [^. The renormalized operator 
0(yu) is defined as the bare operator Obare multiplied by the Z-factor, 

0(yU) = Zoifi] a)Obare{a) . (1) 

For simplicity of presentation we ignore the possibility of operator mixing which would 
require replacing Eq. (|l]) by a matrix equation but have no other effect on the method de- 
veloped in this paper. Here the bare operator Obare(«) is regularized by discretizing it on 
a space-time lattice with lattice spacing a that has units of inverse energy and throughout 
this paper a is assumed to be uniform in all directions. It can also depend on other (bare) 
parameters in the Lagrangian. The renormalized operator 0{n) has low energy matrix ele- 
ments with no dependence on regularization parameters. These low energy matrix elements 
are functions of only physical parameters, such as physical masses, coupling constant, etc., 
defined at the renormalization scale fi, and /i itself. 

The renormalization condition imposed on an operator constructed from a product of n 
quark fields is given by 

Zoifi; a)Z~"/2(/i; a)Toip; a)|p2=^2 = 1 , (2) 
where Zq{fi; a) is the quark field renormalization coefficient, 

qrenifJ') = Zy^ifX] a)qbare{a) ■ (3) 

This specifies our conventions for Zg. The precise prescription for determining Zg is given 
later in Eq. (^. 
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The quantity To{p) is obtained from the amputated matrix element of the operator 
Obare between external off-shell (bare) quark states, in the Landau gauge, by tracing it with 
a projector Pq on a tree-level operator. In the case of quark bilinear operators it is defined 
by the equation 

ro(p; a) = — ^ ^ , (4) 

where 

Ao{p] a) = SbareiV, ay^Go{p] a)SbareiP; Ct)"^ (5) 

is the truncated Green's function. The corresponding equations for the case of four quark 
operators can be found in Section IV D| of this paper. 



The above renormalization condition is independent of the regularization scheme, al- 
though it does depend on a particular choice of the off-shell momenta and the choice of 
gauge. This dependence is not important since the final physical results will not depend on 
this momentum choice, and will be gauge invariant. We always choose the sum of all exter- 
nal momenta to be equal to zero. In the case of the four quark operators all four external 
lines carry momenta that are equal up to an overall sign. 

The renormalization coefficient of the fermion field, Z^^, is defined from the full propa- 
gator. 
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Z,{^;a) = --TriY ^- ] (6) 



and can be calculated numerically using the Ward Identity for the conserved vector current 
V'" in the form 

Kc (p; a) = , (7) 

which implies 

^.(/^;«) = ^Tr(7'^A^c(p;«))^,^^, . (8) 

Substituting this into Eq. for one can see that the RI renormalization condition for 
the conserved vector current is consistent with the equality Zyc = 1. 



III. SCALING TECHNIQUE 

Our scaling technique exploits the renormalizability of lattice QCD which implies that, 
when scaling violations are small, the values of a renormalization coefficient calculated at two 
different values of the lattice spacing will be related to each other by an overall multiplicative 
coefficient independent of the physical conditions such as physical volume or momentum 
scale. 
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Zo{fi,V;a') 
Zo(/i, V; a) 



Roict', a) 



(9) 



We use this property to "sew together" evolutions of renormahzation coefficients over over- 
lapping regions of momentum scale. We are going to describe this technique in detail here. 

In the formula above we explicitly show the dependence of the renormahzation coefficients 
on the physical volume V. While this dependence could be eliminated by taking the infinite 
volume limit, in our scaling method this is not necessary for the intermediate steps where 
the finite volume effects are explicitly taken into account. To avoid complications associated 
with the quark mass renormahzation, the chiral limit is always taken. 

In any Monte-Carlo QCD simulation the size of the lattice in any direction is finite and 
rarely exceeds 32. This puts limits on the momenta accessible in a simulation with one fixed 
value of the coupling constant (3. The components of momenta can assume only a discrete 
set of values. 



For the momenta with large values of Ui the discretization errors become significant. This 
forces one to use momenta with small rii and therefore limits the range over which the 
non-perturbative renormahzation technique outlined above can be carried out. 

Finding a way to study the renormahzation of operators in a wide range of momentum 
scales is important in many applications. One particular but very important example is 
the connection with perturbation theory. Most of the lattice simulations are performed at 
low momentum scales, where the non-perturbative effects play a major role. On the other 
hand, if one needs to relate the results to perturbation theory, for example, to use the 
Wilson coefficients determined from a perturbative calculation, one has to have access to 
large physical momenta (several GeV) where perturbation theory should become accurate. 
This is very difficult to achieve in a simulation at a single value of (3. 

It is natural to think of using the results of numerical simulations at several /3's with over- 
lapping regions of physical momenta with small discretization errors to extend the range over 
which renormahzation coefficients are calculated numerically. The procedure implementing 
this idea is the following (see Figure |ID. One starts at /3(a) which is sufficiently large to in- 
sure that there are momenta well below the lattice cutoff that are in the perturbative regime. 
At such scales the renormahzation coefficients obtained from the Monte-Carlo simulation 
can be reliably related to the perturbative ones. The finite volume effects at this (3{a) can 
be eliminated by numerically taking the infinite volume limit. Alternatively, one can use 
perturbation theory carried out at finite volume. 

The next step involves simulations at a series of /5's with lattice spacing a increasing 
by factors of two. The number of lattice sites in all directions is fixed. The physical 
volume is finite and different in each simulation, increasing by factors of 2^. The range 
of momenta with small discretization errors in each simulation also scales by factors of 
two. The overall range of momentum scales with small discretization errors covered in 
this sequence of simulations can be very large. But the Z-factors obtained with different 
lattice spacings cannot be meaningfully compared directly, since the lattice spacing plays 
the role of the ultraviolet cutoff on the lattice. Changing the cutoff requires complex, 



2n 

Pia = —Hi 



= -(Li-l)/2,...,0,...,Li/2. 



(10) 



5 



generally nonperturbative renormalization of operators in order to keep the physical, low 
energy predictions of the theory unchanged. Therefore, in order to relate the Z-factors 
obtained in simulations at different /?'s, it is necessary to rescale them to be defined with 
the same ultraviolet cutoff equal to the lattice spacing a at some fixed value of (3. To 
achieve this, we use the property Eq. which lies in the heart of our scaling method. For 
the common ultraviolet cutoff we choose, quite arbitrarily, the lattice spacing at /3 = 7.05 
from the first simulation in the series. 

In order to accomplish the program of redefining all of the Z-factors at a common ul- 
traviolet cutoff a it is necessary to know the coefficients Ro{a,a') that relate the Z-factors 
computed at the cutoff a' to the Z-factors that would be found if the cutoff a was used 
instead. To compute the rescaling coefficients Roidya') we perform a series of additional 
simulations. Each additional simulation has the lattice spacing two times larger than the 
corresponding simulation from the main series. At the same time, the number of the lat- 
tice sites in each direction is reduced by a factor of two. With this choice of parameters, 
the physical volume and the values of physical momenta are the same for these two sim- 
ulations. If discretization errors were negligible, the Z-factors from these two simulations 
computed at the same physical momenta would be related by a momentum-independent fac- 
tor i?o(2"'a, 2""*"^a), in accordance with Eq. (^. Here 2"'a and 2"+^a are the lattice spacings 
used correspondingly in the simulation from the main series and the additional simulation, 
while a is the value of the common ultraviolet cutoff chosen to be equal to the lattice spacing 
used in the simulation with the largest (3 in the series. To find this factor we would use 
the ratio Z{fj,,V;2"'a)/Z{fi,V;2"-~^^a) calculated at any value of physical momentum scale 
/i. Here V is the physical volume in both simulations. 

When discretization errors are present, Z-factors, as well as their ratios, have 
momentum-dependent corrections. For the domain wall fermion action these errors are 
quadratic in momentum. 



Zo(/i,\/;2"a) 
Zo(/i,V^;2"+ia) 



i?o(2"a,2"+'a) +0((/i2"a)=') . (11) 



We use quadratic in fia fits to remove the O ((/x2"a)^) term and extract -Ro(2"o, 2"+^a). 
By combining the rescaling coefficients obtained for pairs of lattice spacings, Ro{a,2a), 
Ro{2a,4:a), etc., we can compute the rescaling coefficients Ro{a,2"-a), 

RoicL, 2a) 

Roia, 4a) = i?o(a, 2a)Ro{2a, 4a) (12) 
Ro{a,Sa) = Ro{a,2a)Ro{2a,4:a)Ro{4:a,8a) 



or, more generally, Ro{2^ci,2^a) where n and m are integers. 

With the rescaling coefficients on the left hand side of equations Eq. (|12D we finally can 
apply the same formula Eq. to the Z-factors computed in the main series of simulations 
at different ultraviolet cutoffs to compensate for the non-trivial renormalization required 
when different values of the cutoff are used. In particular, from the original Z-factors in 
the main series, Zo(/i, V^; a), Zo{fJ^,2^V;2a), Zo{fJ^,4:'^V;4:a), etc., we obtain Zo{fJ^,V; a), 
Zoin,2'^V;a), Zo(Ai, 4"^!^; a), etc.. 



6 



Zo{fi,V;a) 
Zo(/i, 2V; a) = Ro{a, 2a)Zo(/i, 2V; 2a) 
Zo(/x, 4V; a) = i?o(a, 4a)Zo(/i, 2V; 4a) 



(13) 



After such a rescaling the ratios of the new Z-factors from all of the simulations are com- 
pletely physical and represent the dependence of the corresponding quantities on the mo- 
mentum scale and volume. We would like to emphasize that direct computation of the 
Z-factors on the left hand side of Eq. (|T3p in the whole momentum scale range for which 
they have been determined would require tremendous computer resources and at the very 
least would be completely impractical. 

After the rescaling, the renormalization coefficients Zq from different simulations evalu- 
ated at the same physical momentum scale will in general differ, because they are calculated 
in different physical volumes. In fact, this method offers one the possibility to study the 
finite volume effects using same size lattices with nevertheless quite different physical vol- 
umes, which can provide certain advantages. The difference caused by the finite volume 
effects gets smaller when the momentum scale becomes large in comparison with the inverse 
size of the volume extent. For some renormalization coefficients these effects are quite large 
(see results for the scalar and pseudoscalar densities and 0^£"^). This effect is illustrated 
schematically in Figure |l| on the graph at the bottom. The labels on the lines show explicit 
dependence on the physical volume, which is different for the simulations at different /?'s. 
If the finite volume effects are small, these lines overlap at the corresponding momentum 
scales. This is (almost) the case for the local and exactly conserved vector and axial cur- 
rents. For any renormalization coefficient there exists a line that corresponds to the infinite 
volume. It is shown as a dotted line on the plot at the bottom of Figure |I|. Since in our 
method we do not need to take this limit, our lines are just approximations of the infinite 
volume line. At each momentum scale the line that has larger physical volume is a closer 
approximation. 

Another reason why the Z-factors computed at the same physical momentum can differ 
between simulations and in fact from their true values is the presence of discretization errors. 
We eliminate these errors from the determination of rescaling coefficients Rq using the 
procedure outlined above, but not from the Z-factors themselves. To remove discretization 
errors from the Z-factors we could use even more additional simulations with matching 
physical volumes so that a continuum limit at each momentum could be taken. The presence 
of these errors does not effect our ability to relate Z-factors between arbitrarily different 
momenta since we only need Rq for that. Since discretization errors are quadratically small 
for small lattice momenta, the values of the computed Z-factors at each physical momentum 
are closer to their true value for data points that come from the simulations at larger (3. 
When the finite volume effects are negligible, so that the expected difference between Z- 
factors at the corresponding physical momenta is small, one can use this observation to 
estimate the magnitude of the discretization error effects. It turns out that their significance 
varies depending on the particular operator for which the Z-factor is computed. As a result 
of our study we found that discretization errors are quite large for the local vector and 
axial currents, and even larger for the exactly conserved vector and axial currents. In the 
16^ volume, for example, only four to five points with the lowest lattice momenta have 
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reasonably small errors. This implies that it is virtually impossible to study the running of 
the renormalization coefficients for these operators without using the scaling technique or 
performing simulations in very large volumes that is extremely expensive computationally. 



IV. NUMERICAL RESULTS 

We use the fermion action which is the hermitian conjugate of the action in Ref. 
where the domain wall fermion method was introduced. The domain wall fermion operator 
that we use is given by 

Dx,s;x',s' = Ss,s'D^lx' + ^x,x'D^,s' (14) 

dIx' = ^ E [(1 - l>:)UxJx+f.,x' + (1 + 7^)f/i',/x-A,x'] + {M, - 4)5,,,, (15) 

[(1 - 75)5.,L.-l5o,s' + (1 + 75)5«,0'5l,-1,s'] (16) 

Operators are constructed from an interpolating operator for the light fermion of the form 

q{x) = Pl*,,o + Pr'^x,l^-i (17) 
q{x) = "^x^oPr + (18) 

Evaluation of the matrix elements is done in Landau gauge. The gauge is fixed on each 
configuration by finding a unitary transformation of the lattice links that maximizes the 
functional 



11=1 



+ UUx)) . (19) 



In our calculations the parameter M5 is set to 1.8. We use lattice volumes 8'^ and 16^ at 
three values of P, 6.0, 6.45, and 7.05. We assume the inverse lattice spacing at /? = 6.0 to be 
1.96GeV, twice that at /? = 6.45, and four times that at (3 = 7.05. These values are based 
on the string tension data in |n|. The scale for (3 = 6.45 was obtained by interpolating 



the data from that paper using the two-loop perturbation formula for the running coupling 
constant. The result for j3 = 7.05 was obtained by extrapolating and is less accurate. 

The chiral limit for the local operators is taken using correlated fits to the results obtained 
at a set of masses mj = 0.004, 0.012, 0.020 in 16^ volume and mj = 0.004, 0.012, 0.020, 0.028 
in 8^ volume. 

In the case of conserved currents all simulations are performed with ruf = 0.05, which is 
small enough for the results to be in the chiral limit. We tested this statement by repeating 
one of the points using mj = 0.02 and found the data to be in agreement within statistical 
errors. 
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A. Local Vector and Axial Currents 



Now we present the results for renormalization coefficients of local vector and axial 
currents. The local vector and axial current densities are defined by 

VT(x) = q^^q{x) (20) 

and 

Al{x) = qYl5q{x) . (21) 

In these equations flavor non-singlet currents are implied. The flavor index is suppressed 
to simplify the notation. The matrix elements of the vector and axial currents in external 
quark states include finite terms that depend on the direction of the external momentum 
in addition to the term proportional to the free-field vertex. In the continuum the form of 
these matrix elements is the following, 

A^4p) = Ay.{p')r + By.ip')^ (22) 

and 

K^ip) = ^^^(/)7"75 + 5^.(^)^75 . (23) 

For the Z-factors we use the definition with the four Lorentz components of matrix elements 
combined by contracting with the gamma matrices, 

^^^^^-TT(AUp;a)r) , , (24) 

and 

lTr(A:;,(p;a)7^5), (25) 



With these definitions the second term in Eqs. (^) and ( ^31) with coefficient -B(p^) is included 



as part of the matrix element. (Alternatively, one could separate the coefficient of the term 
proportional to the free-field vertex, for example, by tracing the /i = component of the 
matrix element calculated at momenta p with p° = with 7° for the vector current and 7*^75 
for the axial current.) 

Figure |^ shows the renormalization coefficients Zq/ZyL and ZqjZj^L versus the energy 
scale in GeV in the chiral limit. The rows correspond to different physical volumes. The 
pairs of points on each plot in the two lower rows are from simulations at pairs of different 
/9's with the lattice spacing differing by factors of two. Since the number of lattice sites 
differs by the inverse of two, the physical volumes are the same. 

Figure ^ demonstrates the procedure of obtaining the rescaling coefficients Ra and Ry 
for two pairs of /3, 6.45, 6.0 and 7.05, 6.45. Each point in Figure |^ is a ratio of the Z-factors 
in Figure ^ determined at the same physical momentum and in the same physical volume. 
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There are four ratios at each pair of /3's, since we only use the lowest four momenta in 8^ 
lattice volume (triangular symbols in Figure 0). 

As described in Section |T|, these ratios would be /z independent in the absence of dis- 
cretization errors. With finite lattice spacing these errors are quadratic in momentum scale 
/i, see Eq. (|11]). In order to remove the effects of discretization errors, we extrapolate the 
ratios linearly in /i^ to yU ^ 0, Figure |^. The rescaling coefficients for the local vector and 
axial density operators turn out to be very close to one. 

Finally, using the rescaling coefficients we rescale the renormalization coefficients accord- 
ing to Eq. (|T3p so that they are all defined with the same ultraviolet cutoff, determined by 
the lattice spacing at /5 = 7.05. The results are shown in Figure ^. There is a pretty good 
agreement between the data computed in different physical volumes at low lattice momenta 
(the first few low momentum points on each graph) that indicates small finite volume effects 
for local vector and axial current data. At the same time, discretization errors start effecting 
data pretty early. To emphasize the actual behavior of the Z-factors we use solid symbols 
for the points that have relatively small discretization errors. These are the same points 
that are used in calculating the rescaling coefficients in Figure 0. 

The renormalization group running of from perturbative analysis [|12| is represented 



by a solid line. This line is defined up to an overall multiplicative constant that we do 
not determine. Instead, we choose it so that the line agrees with our numerical data at 
scales about 10 GeV. The coefficient Zyh should be momentum-independent. It reflects 
the renormalization of the local vector current relative to the exactly conserved on the 
lattice vector current whose renormalization coefficient is equal to one according to the 
Ward Identity, Section [IV C| . The two currents are different at the scales of order a~^, while 
for the scales -C the difference can be absorbed into Zyh. Therefore, the dependence 
of Zql Zyh on the scale [i should be the same as the dependence of Zq itself. There is a good 
agreement between the RG line and the data for Zq/Zyh in the whole range of the scales. 
On the other hand, Zj^l defined by Eq. ( P^D is momentum-independent only at large scales 
11. The reason is the same as for the difference between Zq obtained from the conserved 



vector current and Z^"^) from axial current, see Section |IV C 



B. Scalar and Pseudoscalar operators 

Now we present the results for renormalization coefficients of scalar and pseudoscalar 
operators. The matrix elements of the pseudoscalar density at low energy scales (about 
1-2 GeV) receive contributions from the chiral symmetry breaking in the form of a pole in 
momentum squared and quark mass, 

Zp[p ) p rrif 

The general form of the chiral condensate is given by the equation 

(qq) = a + bmf + (27) 
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The third term in the above equation is due to the fermion zero modes present in our 



quenched simulations, Ref. |T3[- This term, which is inversely proportional to mj, is sup- 
pressed by the square root of the volume and is small in comparison with the constant term 
for the parameters used in our simulations. Therefore, in Eq. (p6D for the pseudoscalar 



(28) 



density the mass dependence of the second term is very well described by the equation 



For the scalar density we have. 



A (r?\ - ^'^^P'^ + ^ ^^^^^ (2^) 



In this case the constant term in the chiral condensate, Eq. (27), is removed by the derivative. 
Therefore, the third term in Eq. (|27| ) which is inversely proportional to rrif, can determine the 
form of the mass-dependence of A5 in small physical volumes. In this case the contribution 
to Asip^) from the second term in Eq. ( PD| ) has the form 

- ■ (30) 

When taking the chiral limit of scalar and pseudoscalar densities at energy scales of a few 
GeV in finite volume, one needs to remove these non-leading, contributions which have 
singular behavior as m/ ^ 0. For a thorough discussion of both of these effects see [0 



The presence of these terms is clearly seen in Figure |] where the data for a 16^ volume 
at /3 = 6.0 for a set of quark masses is presented. Their contribution to the scalar density 
becomes significant only for mj = 0.004 and is almost invisible for mj = 0.02. Figure ^ 
shows results at /3 = 6.45. At this (3 there is little evidence of such small-m/ divergent 
term, which is consistent with the absence of the chiral symmetry breaking in this smaller 
physical volume. 

Figure |^ shows the renormalization coefficients Zg/Zp and Z^/Zs plotted versus the 
energy scale /i in GeV in the chiral limit. To obtain these results, we performed correlated 
fits in m/ at each momentum in the form Eqs. ( ]28|) and (|30D plus finite, mj-independent 
renormalization coefficients Zg/Zp and Zq/Zs- The arrangement of the graphs is the same 
as for the local vector and axial currents. 

The result of the scaling procedure is shown in Figure |^. Solid lines are the renormal- 
ization group running from a 3-loop perturbation theory analysis, [0. They represent the 



evolution of the Z-factors with the energy scale, and are defined up to an overall multi- 
plicative constant which we do not determine. We choose this constant to match the RG 
curve with our data at large scales n around 10 GeV. The data exhibits large finite volume 
effects at small momenta in each of the three different physical volumes. This results in the 
difference between data at the corresponding momenta obtained in different volumes. This 
difference vanishes as the momentum scale becomes large within the individual momentum 
range for each volume. On the other hand, discretization error effects are barely noticeable. 
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C. Conserved Vector and Axial Currents 



The exactly conserved domain wall vector and axial currents can be used to compute Z^. 
Following Ref. [0], with our fermion operator given by Eq. ([I^) these currents are defined 
by the equation 

Jt,{x) = aj J2 J2 J/^{x,s), (31) 

s=0 s=Lj2 

where J = V or A and ay = +1, cta = —1 and 

j,{x, s) = ^,+^,,l±^f/t (a;)vl/,,, - ,l_2:^f/^(a;)vl>,+^,, . (32) 

Here fiavor non-singlet currents are implied with the fiavor indices suppressed for simplicity 
of notation. With the definition of Zg given by Eq. we can use the Ward Identity Eq. (|^ 
for the exactly conserved vector current which implies Eq. (|^). The corresponding Ward 
Identity for the axial current has an additional term that does not vanish in the chiral limit, 
due to the presence of a massless Goldstone boson. It becomes negligible at sufficiently large 
scales /i (see a proof in P] for example). For this reason, when we apply the same equation 
Eq. (^) with A^c, the resulting renormalization coefficient Z^^^ is different from Zg when 
evaluated at scales below ~2 GeV. 

Due to the presence of the fermion fields at separate space-time points in Eqs. (|31| ) and 



(^), propagators must be evaluated for sources at different points on the four-dimensional 
lattice. In addition, the summation over the position in the fifth dimension implies that 
propagators for all of these positions should in principle be calculated. Therefore, the calcu- 
lation of the matrix elements of conserved currents can be very expensive. To minimize the 
amount of computational time, we use a random source estimator to compute the part of 
the sum between s = 1 and s = Lg — 2, with propagators for s = and s = — 1 calculated 
explicitly. Also, instead of calculating all four components of r^(p) for a given momentum 
p, we calculate T^{p) for momenta related to p by permuting its 0th component with the 
rest of the components, 

Tr (7'^^'^(/,p^p^p^)) = Tr {^'T%p' , p\ p' , p')) + Tr {^'V' ip\ p' , p' , p')) 

+Tr (7°^^p^p^p^p^)) +Tr(/^^p^p^p^/)) . (33) 

The computer time required to Fourier transform a matrix element is negligible in compar- 
ison with the calculation of the matrix element itself, therefore the formula above allows us 
to obtain the result with only a quarter of the running time. 

We used mf = 0.05 in our simulations for the conserved vector and axial currents which 
from our experience is sufficiently close to the chiral limit. A test run at /3 = 6.0 with 
rrif = 0.02 with 30 configurations produced results consistent within statistical errors. In 
Figure ^ we plot Zg computed from the conserved vector current using Eqs. (§) and (p3D. The 



same figure shows the renormalization coefficient Zj^^^ obtained by applying these equations 
with the axial current. 
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The data is similar to that obtained for the local vector and axial currents, but shows 
larger discretization errors. In Figure |l^ we plot the evolution of Zq and Z^"^^ factors com- 
puted from the conserved vector and axial currents after rescaling, such that the ultraviolet 
cutoff is equal to the lattice spacing aX (3 = 7.05. The renormalization group running of Zq 
from a perturbative analysis |T2| is represented by a solid line. This line is defined up to an 
overall multiplicative constant that we do not determine. Instead, we choose it so that the 
line agrees with our numerical data at scales about 10 GeV. 

From the results in Figure |10| obtained at different /3's, one can see that while the finite 
volume effects are still small, as in the local current case, the quadratic in deviations due 
to discretization errors that curves the results up become significant very early. Essentially, 
only the first five points with lowest momenta at each (3 seem to have relatively small 
discretization errors. This statement of course is correct only in a 16^ volume. But it does 
suggest that attempts to analyze the momentum dependence of Zq from conserved currents 
using simulations at a single value of /? are subject to large discretization errors. We use 
solid symbols for the data points that have reasonable discretization errors. These are the 
same points that are used in calculating the rescaling coefficients for conserved currents. 

The results from the conserved vector and axial currents agree pretty well at the scales 
of about 2 GeV and above. At the same time, at 2 GeV and below there is large difference 
between Zq and Z^^\ The difference reaches almost 20% at about 1 GeV and is more than 
30% at 0.8 GeV. The numerical results therefore indicate that the additional finite term in 
the axial Ward Identity is large at these scales. The same argument is consistent with the 
data for the rescaling coefficients in Table |. The agreement between rescaling coefficients 
for Zq and Z^^^ is much better for the pair of weaker couplings (3 = 6.45 and P = 7.05 than 
for (3 = 6.0 and f3 = 6.45. 

The results in Figure [1^ are in poor agreement with the renormalization group curve 
from perturbation theory, though it becomes slightly better at large momenta as the scale 
dependence of Zq (and Z^'^'') becomes flatter. This is in sharp contrast with the results for the 
local vector and axial currents. Figure ^ In the absence of scaling violations the results from 
local currents should agree up to an overall momentum-independent factor with the results 
from the conserved currents, as argued in Section [IV A| . We interpret this inconsistency 
between our results for the local and conserved currents as arising from discretization errors 
present in the non-local, conserved currents. These can arise in two ways. First, our 
matching procedure using the small S'^ volumes might not be able to properly account for 
the order 0(a^/i^) terms. This difficulty can be addresses by working on 16^ and 32'^ volumes 
that would allow smaller lattice momenta to be used, therefore reducing the errors. The 
second source of errors may arise from the failure of these non-local conserved currents to 
scale at couplings as strong as (3 = 6.0. 



D. Bk 

In this section we present the results of non-perturbative renormalization of a phe- 
nomenologically important parameter Bk which is relevant to the calculation of the — 
mixing amplitude. The standard definition of Bk is given by 
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... {K'\Ott'mK') _ {K'\Ott'mK') 

This definition implies that Bk renormahzes in the same way as the four-fermion operator 
Ott' = (sr^-^d) (sY^-^d) ^ sr^^ds^ld . (35) 



We closely follow the RI procedure for renormalization of this operators described in |]I6 



In the domain wall fermion formulation the mixing with operators of different chiralities is 
small and is not studied in our paper. For more detailed treatment of this topic the reader 
is referred to Ref. §]. The four-point Green's function is defined by 

Gij,,.(a;i,X2,X3,a;4) = {s{xi)d{x2)0^l='^{Q)s{x'i)d{xi)) . (36) 

For the Fourier transformed Green's function the off-shell external momenta are chosen 
to be equal. Using the upper Roman indices to denote color and the lower Greek indices to 
denote spin, the non-amputated Green's function can be written as 

GBAp)tlT = 2 [{T'{p)t^T'{p)'^f) - (rnp)^f r'^lp)?/ )] , (37) 

where 

r^p)^,^ = 5(p|0)i^7Lp (755(p|0)t75) J . (38) 

5(^10) is a Fourier transformed propagator on a single configuration in Landau gauge. It is 
not translationally invariant. 

The amputated Green's function is calculated from the non-amputated one by multiply- 
ing it with the inverse propagator S{p) averaged over all configurations 

^BAp)ig,T = s-\p)^^:s-\p)^^^'GsAp)i'4Y^^^^ ■ (39) 

The projection operation Eq. (^ is defined by the equation 

1 



(40) 



where Ai'c = 3 is the number of colors. The renormalization condition imposed on the 0^^^^ 
operator is given by the same formula Eq. (^). 

The numerical results for this Z-factor obtained with different values of m/ have no 
statistically significant dependence on ruf. Nevertheless, the same procedure of linear cor- 
related fitting in my was applied to take the chiral limit with the results shown in Figure 
pA| . The final results after rescaling to a new ultraviolet cutoff equal to the lattice spacing at 
j3 = 7.05 are presented in Figure [T^. It is interesting to notice that the finite volume effects 
at low lattice momenta are very significant and for the lowest momenta are of the order of 
50%. This same volume dependence can also be seen from a direct comparison with the 
smaller, 8^ volume at the same p. For illustration purposes we include the result at /3 = 6.0 
in a 8^ volume in Figure O using small circles. The data points for this smaller volume are 
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rescaled using the same factor as the P = 6.0, 16^ volume result, therefore the difference be- 
tween the two results is a direct result of the finite volume effects. The systematic difference 
from the data at /3 = 6.45 that has the same physical volume is due to larger discretization 
errors that are present in the (3 = 6.0 data at the corresponding momenta. 

The solid line on the graph is for the one-loop renormalization group analysis in per- 
turbation theory [|l7l. Since it is defined up to an overall multiplicative constant which 
we do not determine, we matched it with the numerical result at large energy scales of 
about 10 GeV. There is a surprisingly good agreement in the whole range from 1 to 10 
GeV. This might be fortuitous for /i lower than 2 GeV, since at these scales we expect large 
non-perturbative effects. These effects together with the large finite volume effects at low 
momenta could have resulted in the flat shape of numerical results at around 1 GeV. 

V. CONCLUSIONS 

The combination of recursive scaling and NPR offers a promising connection between 
the perturbative and non-perturbative regimes-a connection that is required if perturbative 
errors are to be properly controlled. The results reported here suggest a consistent picture 
for the scaling of renormalization factors over a momentum range between 2 and 10 GeV for 
the case of quenched domain wall fermions for the operators analyzed. The only exception 
to this generally satisfactory picture occurs in the comparison of the local and conserved 
currents, where the larger discretization errors present in the non-local, conserved currents 
may be at fault. To fully resolve this issue and demonstrate the validity of this approach 
more extensive calculations and larger lattice volumes will be required. 
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TABLES 



TABLE L Rcscaling coefficients. 



Z-factor 


R{ai/3 = 6.45), a(/3 = 6.0)) 


i?(a(/3 = 7.05),a(/3 = 6.45)) 


ZqlZyL 


L016(3) 


0.999(2) 


Zq/Z^L 


L014(3) 


0.987(2) 


Zq/Zp 


L088(9) 


1.039(6) 


Zq/Zs 


L046(8) 


1.025(6) 


Zl/Zo^l=^ 


0.881(9) 


0.932(8) 


Zq 


1.188(13) 


1.067(5) 


Zjq 


1.092(13) 


1.047(5) 



TABLE IL Z-factors at /? = 6.0 in 16"^ lattice volume. The mass for the local Z-factors is 
rrif = 0.02, while for the conserved currents (last two columns) rnf = 0.05. 





/i,GeV 


Zq/ ZyL 


Zq/Z^L 


Zq/Zp 


Zq/Zs 


Zl/Zofs=2 


Zq 


^q 


0001 


0.770 


1.1859(91) 


0.9865(71) 


5.610(22) 


1.43(14) 


1.229(25) 


0.914(16) 


0.748(12) 


0011 


1.089 


1.1350(48) 


1.0545(39) 


3.620(13) 


1.430(89) 


1.258(15) 


0.862(13) 


0.790(11) 


0111 


1.333 


1.1142(41) 


1.0731(35) 


2.799(89) 


1.451(65) 


1.265(13) 


0.841(12) 


0.801(12) 


0002 


1.539 


1.0882(27) 


1.0660(26) 


2.361(66) 


1.466(50) 


1.238(10) 


0.824(14) 


0.800(14) 


nil 


1.539 


1.1034(45) 


1.0820(45) 


2.314(54) 


1.456(37) 


1.273(15) 


0.828(11) 


0.802(10) 


0012 


1.721 


1.0915(23) 


1.0749(21) 


2.081(48) 


1.432(36) 


1.249(9) 


0.836(10) 


0.817(10) 


0112 


1.885 


1.0917(25) 


1.0804(23) 


1.907(36) 


1.419(25) 


1.261(9) 


0.841(11) 


0.828(11) 


1112 


2.036 


1.0927(29) 


1.0836(29) 


1.785(28) 


1.405(18) 


1.271(10) 


0.844(12) 


0.836(12) 


0022 


2.177 


1.0913(22) 


1.0835(21) 


1.703(24) 


1.386(16) 


1.268(8) 


0.873(10) 


0.863(10) 


0122 


2.309 


1.0938(24) 


1.0880(24) 


1.639(22) 


1.370(15) 


1.282(9) 


0.888(10) 


0.881(10) 


1122 


2.434 


1.0962(29) 


1.0913(29) 


1.591(18) 


1.364(12) 


1.292(10) 


0.902(11) 


0.897(11) 


0222 


2.666 


1.1004(27) 


1.0966(27) 


1.522(16) 


1.339(12) 


1.307(10) 


0.958(11) 


0.953(11) 


1222 


2.775 


1.1025(32) 


1.0993(32) 


1.497(13) 


1.344(10) 


1.315(11) 


0.977(12) 


0.971(12) 


2222 


3.079 


1.1100(36) 


1.1079(36) 


1.441(10) 


1.328(09) 


1.342(13) 


1.057(14) 


1.053(14) 
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TABLE III. Rescaled ,Z^-factors in the chiral Hmit versus the energy scale. The Z-factors in 



this table arc defined with the ultraviolet cutoff equal to the lattice spacing at f3 = 7.05. The 
physical extent of the lattice in each dimension is 16 times the lattice spacing at /? = 6.0. 



/7,,GcV 




'7 / '7 


Zjq/Zjp 


r/ 1 ry 

Zq/Zs 




Zq 


z\ 


0.770 


1.1910(15) 


1.0045(12) 


2.1989(147) 


2.0848(78) 


1.0188(34) 


1.1591(24) 


0.8554(18) 


1.089 


1.1533(8) 


1.0691(7) 


1.9869(79) 


1.8985(47) 


1.0521(21) 


1.0939(21) 


0.9034(17) 


1.333 


1.1309(7) 


1.0831(6) 


1.8653(55) 


1.7918(37) 


1.0500(19) 


1.0668(20) 


0.9167(18) 


1.539 


1.1014(5) 


1.0726(5) 


1.7911(38) 


1.7010(29) 


1.0202(17) 


1.0447(22) 


0.9153(20) 


1.539 


1.1200(7) 


1.0903(7) 


1.8335(40) 


1.7317(28) 


1.0487(22) 


1.0500(19) 


0.9171(17) 


1.721 


1.1069(5) 


1.0815(5) 


1.7470(31) 


1.6485(23) 


1.0294(16) 


1.0598(18) 


0.9347(17) 


1.885 


1 1076f5'l 


1 0857f5'l 


1 6891 f27') 


1 6035('21') 


1 0389fl6l 


1 0662fl8) 


9466 f 17) 


2.036 


1.1082(5) 


1.0864(5) 


1.6466(24) 


1.5668(21) 


1.0459(17) 


1.0707(20) 


0.9563(18) 


2.177 


1.1074(5) 


1.0880(5) 


1.6230(21) 


1.5363(18) 


1.0445(16) 


1.1067(18) 


0.9870(17) 


2.309 


1.1098(5) 


1.0916(5) 


1.5988(20) 


1.5142(17) 


1.0557(16) 


1.1268(19) 


1.0076(17) 


2.434 


1.1127(5) 


1.0942(5) 


1.5776(20) 


1.4973(17) 


1.0636(17) 


1.1446(20) 


1.0255(18) 


2.666 


1.1169(5) 


1.0989(5) 


1.5500(18) 


1.4673(16) 


1.0755(17) 


1.2148(20) 


1.0895(18) 


2.775 


1.1186(5) 


1.1013(5) 


1.5393(18) 


1.4635(16) 


1.0812(17) 


1.2391(21) 


1.1108(20) 


3.079 


1.1268(6) 


1.1103(6) 


1.5212(18) 


1.4394(16) 


1.1040(18) 


1.3404(23) 


1.2039(22) 



TABLE IV. Rescaled Z-factors in the chiral limit versus the energy scale. The Z-factors in 



this table are defined with the ultraviolet cutoff equal to the lattice spacing at /3 = 7.05. The 



physical 


(wicni of (J 


u' lattice in 


each dimcMision is 16 (Imc 


s lli(^ laldcc s 


")aciug at 3 = 


G.45. 


/i,GeV 


Zq/ZyL 


Zq/Zj^L 


Zq/Zp 


Zq/Zs 




Zq 


^q 


1.539 


1.1052(6) 


1.0861(6) 


2.2889(33) 


2.0970(37) 


1.3771(25) 


0.9897(17) 


0.8713(13) 


2.177 


1.0782(3) 


1.0634(3) 


1.7676(18) 


1.7055(19) 


1.1810(14) 


0.9178(12) 


0.8718(11) 


2.666 


1.0699(3) 


1.0553(3) 


1.5857(14) 


1.5470(15) 


1.1400(13) 


0.8827(12) 


0.8522(11) 


3.079 


1.0595(3) 


1.0457(3) 


1.4952(13) 


1.4609(13) 


1.1045(11) 


0.8665(12) 


0.8419(12) 


3.079 


1.0654(4) 


1.0502(4) 


1.5136(13) 


1.4818(13) 


1.1204(14) 


0.8709(10) 


0.8462(10) 


3.442 


1.0596(3) 


1.0458(3) 


1.4524(11) 


1.4244(11) 


1.1090(11) 


0.8665(9) 


0.8442(9) 


3.771 


1.0609(3) 


1.0478(3) 


1.4101(10) 


1.3840(10) 


1.1142(12) 


0.8667(10) 


0.8461(10) 


4.073 


1.0635(4) 


1.0504(3) 


1.3793(10) 


1.3586(10) 


1.1215(13) 


0.8719(11) 


0.8520(11) 


4.354 


1.0623(3) 


1.0485(3) 


1.3663(9) 


1.3415(10) 


1.1160(11) 


0.8797(9) 


0.8603(9) 


4.618 


1.0654(3) 


1.0517(3) 


1.3484(9) 


1.3261(9) 


1.1264(12) 


0.8906(10) 


0.8717(9) 


4.868 


1.0684(4) 


1.0550(4) 


1.3336(9) 


1.3135(9) 


1.1358(13) 


0.9050(11) 


0.8862(11) 


5.333 


1.0717(3) 


1.0580(3) 


1.3150(9) 


1.2969(9) 


1.1466(13) 


0.9365(10) 


0.9176(10) 


5.550 


1.0751(4) 


1.0613(4) 


1.3110(9) 


1.2904(9) 


1.1575(14) 


0.9572(12) 


0.9382(12) 


6.158 


1.0823(4) 


1.0688(4) 


1.2990(8) 


1.2806(9) 


1.1798(15) 


1.0215(13) 


1.0018(13) 
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TABLE V. Rescaled Z-factors in the chiral limit versus the energy scale. The Z-factors in this 



tabic arc defined with the ultraviolet CTitoff equal to the lattice spacing at /? = 7.05. The physical 
extent of the lattice in each dimension is 16 times the lattice spacing at /3 = 7.05. 







^ql 




^q/^S 


'72 17 

'^ql^Off^'^ 


Zq 




3.079 


1.0748(8) 


1.0624(4) 


1.8674(28) 


1.8050(30) 


1.5200(28) 


0.8407(12) 


0.8100(12) 


4.354 


1.0630(2) 


1.0604(2) 


1.4970(13) 


1.4777(13) 


1.2637(12) 


0.8324(8) 


0.8247(8) 


5.333 


1.0593(2) 


1.0586(2) 


1.3698(9) 


1.3612(9) 


1.2041(10) 


0.8074(8) 


0.8032(8) 


6.158 


1.0416(1) 


1.0406(1) 


1.3058(8) 


1.3017(8) 


1.1472(5) 


0.7888(7) 


0.7848(7) 


6.158 


1.0571(3) 


1.0572(3) 


1.3208(7) 


1.3176(7) 


1.1850(11) 


0.7958(8) 


0.7931(8) 


6.884 


1.0452(1) 


1.0442(1) 


1.2806(5) 


1.2752(5) 


1.1504(5) 


0.7883(7) 


0.7860(7) 


7.541 


1 0494f2') 


1 0487('2) 


1 2527f5') 


1 2492 ('5) 


1 1589(5) 


7854('6') 


7837f6) 


8.146 


1.0518(2) 


1.0517(2) 


1.2320(5) 


1.2291(5) 


1.1649(7) 


0.7877(7) 


0.7862(7) 


8.708 


1.0483(1) 


1.0475(1) 


1.2238(4) 


1.2199(4) 


1.1546(4) 


0.7885(6) 


0.7872(6) 


9.236 


1.0525(2) 


1.0519(2) 


1.2127(4) 


1.2100(4) 


1.1663(6) 


0.7954(6) 


0.7942(6) 


9.736 


1.0562(2) 


1.0563(2) 


1.2033(4) 


1.2012(4) 


1.1794(8) 


0.8053(7) 


0.8043(7) 


10.665 


1.0590(2) 


1.0587(2) 


1.1941(4) 


1.1925(4) 


1.1861(7) 


0.8247(6) 


0.8239(6) 


11.101 


1.0622(3) 


1.0621(3) 


1.1896(4) 


1.1882(4) 


1.1969(9) 


0.8400(7) 


0.8392(7) 


12.315 


1.0678(3) 


1.0676(3) 


1.1806(4) 


1.1790(4) 


1.2139(11) 


0.8831(8) 


0.8824(8) 
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FIGURES 



ZJii,4 V;4a) 




(Vi, V;a) 



|i,GeV 



Zq(|li, V; a) 

R (a,2a) = Urn r^^^ ^ \ 



16' 



|i,GeV 



Z^(^ , 4 V; a) = RJaM) ZJ \i, 4 V; 4a) 



(^i , r V; a) = R (fl,2fl) (|li , 2^ 2a) 




Z^(|I,oo; a) 



|Li,GeV 



FIG. 1. Scaling procedure. Panel 1 shows Z-factors obtained at a series of different /3's and 
same size lattices corresponding to different physical volumes. These Z-factors cannot be mean- 
ingfully compared directly since they are defined with different values of lattice cutoff. Panel 2 
shows the procedure of obtaining the rescaling coefficients Rq- Panel 3 shows Z-factors rescaled 
using Rq and therefore defined with the same lattice cutoff a. The ratios of the Z-factors on this 
graph have direct physical meaning and arc cutoff-independent. Since different physical volumes 
are used, these Z-factors do not necessarily overlap at the corresponding momenta. 
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FIG. 2. Zq/ZyL and Zq/ZyL in the chiral limit from local vector and axial currents. Each pair 

of graphs is for the same physical volume. 
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FIG. 3. The ratios of the renormahzation factors Zq/ZyL and ZqjZj^L obtained at different /?'s 
in the same physical volumes versus momentum squared. The deviation from constant values are 
due to order (/xa)^ discretization errors. Linear fits in [i} are used to find the rescaling coefficients 
between renormalization factors at different /3's. 
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FIG. 4. Zq/ZyL and Z^/Zj^l for local vector and axial currents in a wide range of scales /U. All 

graphs have the same ultraviolet cutoff equal to the lattice spacing at (3 = 7.05. Different symbols 
correspond to data in different physical volumes (labeled by /? at which they are obtained). Solid 
symbols are for data with small discretization errors. 
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FIG. 5. Mass poles in the pseudoscalar and scalar operators in a 16^ lattice volume at /3 

These poles have to be subtracted at each momentum. 
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FIG. 6. Mass pole effects are mucfi smaller at /3 = 6.45 in 16^ lattice volume. 
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FIG. 7. Zq/Zp and ZqjZg from pseudoscalar and scalar density after subtraction of mj poles. 
The graphs from top to bottom correspond to different physical volumes, while each pair of graphs 
is for the same physical volume. 
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FIG. 8. Zq/Zp and Zq/Zs in a wide range of scales /U. All graphs have the same ultraviolet 

cutoff equal to the lattice spacing at /? = 7.05. Different symbols correspond to different physical 
volumes (labeled by the P at which they are obtained). The finite volume effects are noticeable for 
momenta of the order of the inverse box size. 
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FIG. 10. Zg and plotted for a wide range of scales fj,. All graphs have the same ultraviolet 

cutoff equal to the lattice spacing at /3 = 7.05. Different symbols correspond to data in different 
physical volumes (labeled by the /3 at which they are obtained). Solid symbols are for data with 
small discretization errors. 
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FIG. 11. Zg/ZQAs=2 in the chiral limit from local vector and axial currents. The graphs from 
top to bottom correspond to different physical volumes, while each pair of graphs is for the same 
physical volume. 



30 



1.6 



<> 



C) 

[] 



Op=6.0 
□ p=6.45 
Op=7.05 

o |3=6.0, 8^ volume 



O 




|i, GeV 

1 10 

FIG. 12. Z'^/Zqas=2 plotted for a wide range of scales ji. All graphs have the same ultraviolet 
cutoff equal to the lattice spacing at (3 = 7.05. Different symbols correspond to data in different 
physical volumes (labeled by the f3 at which they are obtained). The finite volume effects are 
large for momenta on the order of the inverse box size. Large symbols correspond to the 16^ 
lattice volume. Small circles denote data obtained in a 8*^ volume at /3 = 6.0. Since no relative 
renormalization is involved for the two calculations at (3 = 6.0, the difference between the points 
is due to the finite volume effects only. The difference between the small circles and squares is 
due to discretization errors at the corresponding momenta. These errors are larger for the results 
obtained at the smaller /3 = 6.0. 
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